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Abstract 

The PIENU experiment at TRIUMF is aimed at a measurement of the branching ratio 

R e /d~Y^(ji + —> e + v e ) + ( 7r + — > e + V(.y))/r((7r + — > /r + v^) + (n + —> ^v^y)^ with precision <0.1%. Incident pions, deliv¬ 
ered at the rate of 60 kHz with momentum 75 MeV/c, were degraded and stopped in a plastic scintillator target. Pions 
and their decay product positrons were detected with plastic scintillators and tracked with multiwire proportional 
chambers and silicon strip detectors. The energies of the positrons were measured in a spectrometer consisting of a 
large Nal(Tf) crystal surrounded by an array of pure Csl crystals. This paper provides a description of the PIENU 
experimental apparatus and its performance in pursuit of R' d ‘. 
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1. Introduction 

Measurement of the branching ratio R ellJ 

r((;r + -> e + v e ) + (7r + -> e + v e yj) 

r((7T + -4 jU+Vp) + ( n + -> /U+V^y)) 

compared to its precise Standard Model (SM) prediction provides a stringent test of the hypothesis of e-/j universality 
in weak interactions. The most recent experimental measurements of the branching ratio are 

A™ = (1-2265 ± 0.0034(sfaf) ± 0.0044(iysf)) x 10" 4 QM 
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and 


R*Jf = (1.2346 ± 0.0035 (stat) ± 0.0036 (syst)) x 10" 4 0. 

These measurements are in agreement with the prediction: 

R& = 1.2352(1) xHT 4 00. 

The order of magnitude gap in the precision of the experimental measurements and the theoretical prediction provides 
an opportunity to test the SM and to search for physics beyond it. Because of helicity-suppression of 7r + —> e + v e decay 
in the SM, helicity-unsuppressed contributions due to pseudoscalar or scalar couplings not present in the SM could 
result in a deviation of R'^‘ from the SM prediction. A measurement of R en ‘ with <0.1% precision can probe new 
physics at mass scales up to 1000 TeV for pseudoscalar interactions 0. 

The PIENU experimental technique is based on observing positrons from decays of pions at rest in an active 
plastic scintillator target. Muons from the dominant n + —> fJ 1 vy, decay, which have kinetic energy of 4.1 MeV and 
range of about 1 mm, were stopped in the target. By measuring the time and energy of positrons from the decays 
7T + —> e + v e and /.i + —» e + v e v ll following the decay n + —» /r + v A , (n + —» /r + —> e + decay chain), R e/fl can be obtained 
after applying several corrections. Since the positrons from the 7r + —» e + v e decay are monochromatic {T n ^„, y = 69.3 
MeV), and positrons from the n + — * /i + —> e + decay chain have a three-body decay energy spectrum with highest 
positron kinetic energy cutoff at 52.3 MeV, the measured energy is used to separate the two decay modes. The raw 
branching ratio is extracted by performing a simultaneous fit of the time spectra of high energy (T e >52 MeV) and 
low energy (T e <52 MeV) samples to determine the yields of n + —> e + v e and 7r + —» /r + —» e + decay chain events 
while taking backgrounds into account. Normalization factors, such as the solid angle of positron detection, cancel to 
first order, and the measured decay yields must be corrected for small energy-dependent effects, such as those from 
multiple Coulomb scattering, Bhabha scattering, and positron annihilation. For the final branching ratio determination 
additional corrections are applied including one which accounts for those n + —» e + v e decays excluded from the high 
energy time spectrum due to electromagnetic shower leakage and other effects. The latter correction and its systematic 
uncertainty were empirically determined from data and measurements of the response of the calorimetry system. 

The PIENU detector had a solid angle acceptance^for positrons emitted from the stopping target of 25%, which is 
eight times higher than the previous TRIUMF experiment 0. In addition to measuring R e ^, because of the increased 
acceptance and better energy resolution, the PIENU experiment was also able to perform a more sensitive search for 
heavy neutrinos in the n + —> e + vn decay 0. 

2. Experimental overview 

The TRIUMF cyclotron delivered a 500 MeV proton beam with an intensity of 120 /rA to a 12 mm thick beryllium 
pion production target. A positively charged beam of momentum P = 15 MeV/c was collected at 135° and transported 
in vacuum in the upgraded M13 channcj^][[8) to the PIENU apparatus. 

The PIENU apparatus, schematically shown in Figure |TJ was placed immediately downstream of the exit of the 
M13 beamline. The end of the vacuum pipe was covered with a 101.6 mm diameter, 72.6 /rm-thick stainless steel 
window. To reduce the neutral background coming from the pion production target area and collimators in the beam¬ 
line, a 24.8 cm-thick 2.57 x 2.29 m 2 (width x height) steel wall with a hole (26.67 cm diameter) for the beam pipe to 
go through was installed upstream of the detector assembly. 

During normal data taking, the beam was composed of 84% n + , 14% /i + and 2% e + . It was imaged by a circular 
multiwire proportional chamber package (WC1 and WC2, containing three wire planes each). Following WC2, the 
beam was degraded by two plastic scintillator counters B1 and B2 (beam counters) used for measuring time and 
energy loss for particle identification. The beam counters were followed by two pairs of Si strip detectors (SI and S2) 
with strips oriented along the X and Y axe^j An 8 mm thick plastic scintillator target (B3) was sandwiched between 


5 Acceptance was defined by two plastic scintillator counters located downstream of the pion stopping target. 

6 The M13 beamline at TRIUMF was modified in order to reduce the positron contamination in the beam. 

7 A right-handed coordinate system is used with the origin at the center of the target counter. The Z-axis points in the beam direction, and the 
Y-axis points upwards. 
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Table 1: PIENU detector parameters. 


Plastic scintillator counters 

Trigger Counters 

B1 

B2 

B3 

T1 

T2 

Size in X (Inner radius) 

100 mm 

45 mm 

70 mm 

80 mm 

(0) mm 

Size in Y (Outer radius) 

100 mm 

45 mm 

70 mm 

80 mm 

(171.45) mm 

Size in Z 

6.604 mm 

3.07 mm 

8.05 mm 

3.04 mm 

6.6 mm 

Z position 

-39.03 mm 

-30.02 mm 

0 mm 

19.92 mm 

72.18 mm 

Photomultiplier model/ 

H3178-51 

83112-511 

XP2262B 

83112-511 

H3165-10 

manufacturer 

Hamamatsu 

Burle 

Photonis 

Burle 

Hamamatsu 

Photocathode diameter 

34 mm 

22 mm 

44 mm 

22 mm 

10 mm 

Veto Counters 



VI 

V2 

V3 

Inner radius 



40 mm 

107.95 mm 

177.8 mm 

Outer radius 



52 mm 

150.65 mm 

241.3 mm 

Size in Z 



3.175 mm 

6.35 mm 

6.35 mm 

Photomultiplier model/ 



H3164-10 

H3165-10 

Photomultiplier manufacturer 


Hamamatsu 

Hamamatsu 

Photomultiplier Photocathode diameter 


8 mm 

10 mm 

Tracking detectors 

Multi Wire Proportional Chambers 


WCl 

WC2 

WC3 

Wire spacing 



0.8 

mm 

2.4 mm 

Number of Planes/wires/readout channels 


3/120/40 

3/96/48 

Active area diameter 



96.0 mm 

230.4 mm 

Cathode to Anode spacing 



1.6 

mm 

2.0 mm 

Wire diameter 




15 fim 


Wire orientation 




0°,+120°,-120° 


Silicon Strip Detector Pair (X and Y oriented strips) 



S1/S2/S3 

Active area 





61 x 61 mm 2 

Silicon strip pitch 





80 fim 

Effective pitch after binding 4 strips 




320 fim 

Number of planes/readout channels per plane 



2/48 

Thickness (size in Z) 





0.285 mm 

Separation between X and Y strip detectors 




12 mm 

Electromagnetic calorimeter 

Crystal 




NaI(T0 

Csl 

Number used 




1 

97 

Energy Resolution (FWHM) at 70 MeV 



2.2% 

-10% 

Thickness (size in Z) 




480 mm 

250 mm 

Outer Radius 




480 mm 

- 

Approximate width x height for pentagon shaped Csl crystals 

- 

90x 80 mm 2 

Number of PMTs per crystal 



19 

1 

Hamamatsu PMT model (central PMT for Nal(Tf) was 

R1911-07) 

R1911 

R5543 

Photomultiplier Photocathde diameter 



76.2 

mm 


3 



Y 


Csl Rings 


V3 


10 cm 
i-1 


WC3 


V 

V2 [ 


WC1/2 


VI 


Beam 


B1 


B2 


B3 


T1 


VI 



SI S2 S3 


V2 [ 
V3 



Figure 1: Schematic diagram of the experimental apparatus. See text. 


S2 and another X-Y Si strip pair (S3). A circular multiwire proportional chamber WC3, containing three wire planes, 
was positioned downstream of S3 and sandwiched between two plastic scintillators, T1 and T2 (telescope counters). 

Incoming pions stopped in the center (+1 mm along the Z-axis) of B3, and then decayed at rest. A coincidence 
of signals from T1 and T2 was used to define the on-line positron acceptance. WC3 and S3 were used to reconstruct 
charged particle tracks and define the final acceptance off-line. The 15% solid-angle acceptance of the crystal spec¬ 
trometer was defined by a radius cut (r=6 cm) in the central plane of WC3, which was located 5.5 cm downstrearr ^of 
B3. Accepted decay products entered a large single crystal Nal(Tf) detector (48 cm long and 48 cm diameter) (9) . In 
order to absorb shower leakage from the Nal(Tf) crystal and to detect large angle photons from radiative n + —> e + v e 
y decay, it was enveloped by four rings of pure Csl crystals M1 Oi l. The Nal(Tf) and Csl crystal array served as the 
calorimeter to measure the energy of the decay product positrons. 

Three annular scintillators, which covered inactive material, were used for systematic studies: VI covered the 
front frame of WC1, and V2 and V3 covered the frame of WC3 and the front flange of the NaT(Tf) crystal enclosure, 
respectively. Detector parameters are given in Table [T] 


3. Design and construction 

The experimental apparatus was constmcted in modules (PIENU-I and PIENU-II) for simplification of installation 
and servicing. PIENU-I consisted of VI, WC1, WC2, Bl, B2, SI, S2, B3, S3 and Tl, while PIENU-II included WC3, 
T2, V2, V3, Nal(Tf) and Csl. 


8 Distance from the center of B3 to the central plane of WC3 
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Figure 2: Computer generated rendering of the PIENU detector and the steel wall used as passive shielding at the end of the M13 beamline. 


PIENU-I was attached to the end of the beam pipe and could be divided into three sub-assemblies (V1+WC1+WC2, 
B1+B2+S1+S2 and B3+S3+T1). The PIENU-II detector assembly was supported on a cart. During normal running 
conditions the cart rested on rails which were aligned with the beamline and PIENU-I in the X-Y plane. Modularity 
allowed a number of supplemental measurements to be performed aimed at evaluating systematic effects during data 
taking. One of the measurements was determination of the response of the Nal(Tf) and Csl crystal array to a beam of 
70 MeV/c incident positrons at various positions and angles of entrance as described in Section [573| This measurement 
was made possible by leaving only the first sub-assembly of PIENU-I in place, containing V1+WC1+WC2, which 
detected the incoming beam, and removing the rails in order to allow the PIENU-II assembly to be moved freely 
within the experimental area. 

In 2010, the experimental area was enclosed in a temperature-controlled tent in order to make temperature-related 
changes in detector’s performance negligible. The temperature inside the tent was kept at 20° ± 0.5° C. 

4. Scintillation counters 

The plastic scintillators were fabricated from Bicron BC-408 (polyvinyltoluene) scintillator lUTI chosen for its 
high light-output, fast rise (0.9 ns), and short decay (2.1 ns) times. The B1 beam scintillator was large enough to fully 
cover the aperture of the multiwire chamber package WC1 and WC2. The B2 scintillator, which was smaller than the 
B3 scintillator, ensured that the particles going through B2 would hit B.f] Tl, located downstream of B3, defined the 
positron timing. Due to space restrictions and their circular shape, T2 and the veto scintillators (V1-V3) were read 
out with 1 mm diameter wave-length-shifting (WLS) fibers (Kuraray Y-ll) embedded in grooves (1.lx 1.1 mm 2 in 
cross-section) machined on the scintillator surfaces; groove spacing was 10 mm, 3 mm, 10 mm and 11 mm for T2, 
VI, V2 and V^jcounters respectively. 


9 The rectangular-shaped B3 and Tl counters were rotated by 45° around the beam axis due to spatial limitations of the mechanical assembly 
which supported the scintillator along with their respective light guides, photo multiplier tubes (PMTs), and a number of other subdetectors (WC1, 
WC2, SI, S2, S3). 

10 T2 had linear parallel grooves traversing the counter’s surface, while annular veto counters VI,V2 and V3 had circular grooves. 
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Figure 4: Schematic drawing of the T2 scintillator. 33 parallel grooves 
(1.1 x 1.1 ram 2 ) machined on the surface of the counter were made 
deeper towards the edge to allow the WLS fibers to exit perpendicular 
to the counter’s surface (see text for details). One such readout WLS 
fiber is shown for illustration. 


Figure 3: Schematic view of the B1, B2, B3 and T1 scintillator readout. 
Light guides were bent adiabatically and coupled to PMTs through 
acrylic cylinders. 


Collection of the scintillation light for all the rectangular-shaped scintillators was done with adiabatic light guides 
which were cut from UV transparent acrylic sheets and polished. It was important to have a uniform light collection 
efficiency in order to maximize energy resolution. A schematic of the four-way readout configuration is shown in 
Figure [3] Simulations of the optical transport with Detect2000 fl2l were used to predict the uniformity of the light 
collection in this configuration, and were later confirmed with measurements. For example, a variation of < 1.2% was 
observed in the visible energy deposited by beam muons in B1 over the counter’s surface. 




2 2.5 3 3.5 4 4.5 5 

Target Energy [MeV] 


Figure 5: The observed time difference between T2 and Tl. The T1 
time is determined by fitting a waveform, while the T2 time is defined 
as the time of the earliest pulse found within the waveform. Due to 
WLS fiber readout of T2 several pulses could be seen for one particle. 
Histogram is data, and the smooth line is a Gaussian fit. 


Figure 6: Visible energy deposited in B 3 by a /j + with kinetic energy 
of 4.1 MeV from n + —» fi + v^ decay. The data is represented by points, 
and a Monte Carlo simulation by the histogram. The mean visible en¬ 
ergy is approximately 3.2 MeV due to saturation of scintillation light. 
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Each scintillator was read out with four photo-multiplier tubes (PMTs), except for VI, V2 and V3 which were 
each read out with a single PMT. PMTs were selected by photocathode size and ability to match either the 425 nm 
wavelength of the maximum emission of the Bicron-408 or the 476 nm wavelength of the emission peak of the WLS 
fibers. Information on the scintillator size and PMTs used is listed in Table[l] All plastic scintillator PMTs were read 
out with 500 MHz ADCs (see Section fO} . 

T2 was the largest scintillator read out with WLS fibers, having 33 parallel grooves machined on the downstream 
side. Near the edge of the counter, the grooves were made gradually deeper so that the WLS fibers lying in the grooves 
could exit perpendicular to the counter’s surface with minimum transition space and point upstream where they could 
be bundled and read out by PMTs (see Figure[4]for schematic of the T2 grooves). The minimum radius of the groove 
curvature at the edges of T2 was 19.5 mm. Libers were thermally pre-bent to reduce stress. T2 fibers were read out on 
both ends, and bundled into two groups (1-16 and 17-33). As a result, T2 was read out using four PMTs. It had a time 
resolution of <x t =1.4 ns for minimum ionizing particles when measured against T1 (see Figure [5]> and its efficiency 
for positrons was measured to be greater than 99.7% within the acceptance region. 

The most critical scintillation counter, B3, had an energy resolution of 5% r.m.s. ( T V i S jbi e = 3.22 + 0.16 MeV) for 
4.1 MeV fi' from n + —> /j ' decay (see Figure [6jl. The time difference between T1 and B1 was used to determine 
the decay time with a resolution of cr, =138 ps for beam positrons. Both T1 and B1 times were determined from fits 
of 500 MHz-sampling-frequency waveforms. 


5. Calorimeter 

The purpose of the PIENU calorimeter system was to detect decay positrons and photons from the 7r + —» e + v e and 
7T + —» e + v e y decays as well as from the /u + —> e + v e v fJ and /z + —> e + v e v li y following the decay n + —> It was 

composed of a large monolithic Nal(Tf) crystal surrounded by 97 pure Csl crystals. Scintillation light collection was 
realized through directly coupled PMTs. 

5.1. Nal(Tl) crystal 

The main calorimeter element was a single crystal of thallium-doped Nal on-loan from Brookhaven National 
Laboratory |9). The crystal was enclosed in a 3 mm thick aluminum enclosure which had nineteen 76.2 mm diameter 
circular quartz windows at the rear end. The aluminum front face of the Nal(Tf) was 0.5 mm thick. Each window 
was viewed by a 76.2 mm diameter Hamamatsu R1911 PMT (with the exception of the central PMT R1911-07) 
surrounded by a /./-metal shield. 

The surface of the crystal was covered with reflective material. An optical simulation was performed using De- 
tect2000 UH to study the dependence of the energy deposited in the crystal on the entrance location of the particle. 
Light emitted by the crystal was uniformly rcflccteq 11 [ and thus a similar amount of light was seen by individual 
PMTs independent of the entrance position of the ionizing particle on the front face of the Nal(Tf) crystal. This 
was confirmed within 2% by bench tests using 511 KeV and 1274 KeV y rays from a 22 Na radioactive source. The 
calibration of the Nal(Tf ) crystal was done with 75 MeV/c beam positrons using a dedicated trigger. 

5.2. Csl crystal array 

Each crystal within the Csl array was 25 cm (13.5 radiation lengths) long and had pentagonal shape, with average 
width of 9 cm and height of 8 cm (4.5 radiation lengths). The crystals were arranged to form two upstream and two 
downstream concentric layers around the Nal(Tf). Each concentric layer of crystals was supported by a 2 mm thick 
stainless steel cylinders with a 2-mm thick fin separating and supporting every 3-5 crystals from its neighbors (6 fins 
per layer). The resulting array of Csl crystals was 50 cm long and 16 cm (9 radiation lengths) thick in the radial 
direction. The Csl array was continuously flushed with nitrogen gas to keep the humidity level low. Each crystal 


"For the simulation both specular and diffuse reflectors were modeled, with the reflection coefficient varying between 0.95 to 0.99 and a 
refractive index of 1.85 was assumed for the Nal('IY’) crystal. 
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was read out by a fine-mesh 76.2 mm diameter Hamamatsu R5543 PMT lITS p^ Thc pure Csl pulse shape has two 
components: fast and slow with 30 ns and 680 ns decay times, respectively. The fast component constitutes about 
20% of the total pulse. In order to suppress the slow component, a UV-transmitting optical filter was present in front of 
each PMT. Each crystal had a Yal 03 :Ce 245 Am source ED attached to its front face to monitor the crystal’s light output 
and the PMT’s gain. The source emits light at a frequency of 50 Hz with similar wavelength and pulse width as the Csl 
scintillation light with an equivalent energy deposit of 10 MeV. Each crystal was also connected via a quartz fiber to 
the output of a Xe lamp whose flash was triggered at 2 Hz during data taking. This Xe lamp system was implemented 
to monitor Csl PMT’s gains. Seven reference PMTs of the same type as the ones used for the calorimeter were 
enclosed in an incubator maintained at a constant temperature of 24.0 °C. The Xe lamp was enclosed in an identical 
incubator at the same temperature. The Xe-lamp was connected to the reference PMTs through the same system as 
the Csl crystals so that any changes in the Xe-lamp output could be tracked by the reference PMTf 2 ] 

The energy calibration of the Csl array was performed with cosmic rays. A dedicated trigger, based on coincidence 
of the two concentric Csl rings, provided a sample of cosmic ray data. The energy peak due to minimum ionizing 
muons going through the crystals was compared with a detailed Monte Carlo (MC) simulation of the system to obtain 
the energy calibration. 




Figure 7: Left: Response of the Nal(TY) crystal to a positron beam of momentum 70 MeV/c. Below the full energy peak, the additional low energy 
peak structures are due to photoabsorption followed by neutron escape from the crystal p6) . Right: n + —» e + v e decay spectrum with additional 
cuts for suppressing the n + —> n + v^ decays, which can be seen in the energy spectrum up to approximately 50 MeV. The first structure due to 
photoabsorption in the crystal becomes visible to the left of the main peak (from the n + —» e + v e decay positron at around 65 MeV) as one restricts 
the acceptance. The shaded histogram represents a tighter-than-nominal 40 mm radial acceptance cut and the unfilled histogram represents the 
nominal acceptance cut of 60 mm. 


5.3. Performance 

Nal(TY): the response of the Nal(TT) crystal was investigated by reconfiguring the detector to include PIENU-II 
and only the first sub-assembly of PIENU-I (VI, WC1, and WC2). 

The response of the crystal to a 70 MeV/c positron beam entering its front face at the center is shown in Figure [7] 
(left). The high energy peak corresponds to the full energy of the beam positrons. After deconvoluting the beam mo¬ 
mentum width of 0.5%, the resolution at 70 MeV/c was 2.2% (FWHM). Below the beam energy, additional structures 


12 The components of the magnetic background in the experimental area have been measured to be <2 Gauss at the location of the detector, well 
within the operational specifications of these fine-meshed PMTs. Before being brought to TRIUMF, the crystals and their PMTs were used as 
endcap photon-veto detector in the E949 experiment at BNL E). 

13 The Xe lamp system was found to be stable to less than 1% level and was used for monitoring. Since the whole experimental apparatus was 
enclosed in the temperature controlled tent, no temperature-dependent variations in the PMT gains were observed, while the PMT gain degradation 
over time was on the order of «1 % per year. 
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are visible flT6ll due to nuclear photoabsorption and subsequent neutron escape. In this process, one or more photons in 
the shower initiated by the positron are absorbed by iodine and neutrons emitted after the absorption escape the crystal 
reducing the measured energy by the neutron binding energy of approximately 8 MeV. Structures corresponding to up 
to three escaping neutrons can be seen in Figure [7] (left). This mechanism was consistent with a MC simulation of the 
crystal response including electromagnetic and hadronic interactions Da. 



Figure 8: Fraction of events below 50 MeV for the Nal(TY) crystal alone and the Nal(TT) plus Csl combination as a function of the angle between 
the beam and the calorimeter horizontal axis . Error bars are too small to be visible. 


Photonuclear effects are also present when measuring the n + —» e + v e decay positrons. An energy spectrum of the 
7r + —* e + v e decay positrons with a small (<15%) contamination of positrons from the n —* /r —» e decay chains is 
shown in Figure [7] (right), where the presence of a structure due to one-neutron escape below the main peak (see also 
0) is evident. 

Csl Array: The purpose of the Csl array was the collection of the electromagnetic shower leakage from the 
main Nal(TC) calorimeter. A critical part of the PIENU experiment was the reduction and measurement of the low 
energy tail from the n ' —> e + v e decay to reduce systematic uncertainties in the correction for events buried under the 
n —> n —> e decay chain energy distribution. During the measurement of the calorimeter response a positron beam 
was used to investigate the fraction of the low energy tail as a function of the angle between the beam and the Nal(Tf) 
crystal face. The crystal face was rotated to 10 different angles (up to 48°) with respect to the beam. The results are 
shown in Figure [8] where the fraction of events below 50 MeV (tail fraction) is reported as a function of the angle. 
When only the energy of the Nal(Tf) crystal is considered, the low energy tail fraction reaches 12% for the highest 
angle, while if the Csl energy is added, the fraction never exceeds 2% over the whole angular range. 

6. Tracking devices 

The beam wire chambers WC1 and WC2 were used to track the incoming pions (see Figure|9|. Together with the 
two X-Y planes of silicon strip detectors (SI and S2) they could detect pion decay-in-flight events upstream of the 
target (see jj |6.3| l. The downstream trackers (S3 and WC3) were used to track the decay positrons after the target. The 
reconstructed track’s position at the center plane of WC3 was used to define the acceptance of the positrons at the 
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entrance of the calorimeters. Upstream (WC1, WC2, SI, and S2) and downstream (S3 and WC3) tracker information 
was also used to reconstruct the position of the pion’s decay vertex in the target. WC3 was also used to detect pile-up 
events which were not detected by the upstream trackers or S3 which had smaller recording time windows (see Section 

im 


6.1. Multiwire proportional chambers 

Each wire chamber (WC1, WC2 and WC3) consisted of three wire planes which were rotated by an angle of 120° 
with respect to each other. The chambers were filled with a tetrafluoromethane (CF 4 ) - isobutane (C 4 H 10 ) mixture 
(80% - 20%) at atmospheric pressure. Each wire plane of WC1 and WC2 had 120 wires (0.8 mm in pitch) grouped 
by three. Each group was connected to a read-out pad. The number of read-out channels was 40 per plane. The 
active diameters of WC1 and WC2 were 9.6 cm. The chamber WC3, mounted on the flange of the Nal(Tf) crystal 
enclosure, measured the position of the decay positrons at the entrance of the calorimeter, and defined the acceptance 
region. Each plane had 96 wires with a pitch of 2.4 mm. The wires were grouped in pairs which reduced the number 
of channels to 48. The active diameter of WC3 was 23 cm. The efficiency of every plane was measured to be > 99% 
for beam positrons. 



Figure 9: Beam profile plots. Scatter plot (Y vs X) of the beam profile at the center of the WC1, the red box indicates the typical boundaries of 
the incoming particle selection cut (left plot). X and Y profiles of the beam at the center of WC1 (center and right plots, respectively). X and Y 
coordinates are determined by reconstructing the incoming particle’s track using WC1 and WC2. 


6.2. Silicon strip detectors 

Each set of silicon strip detectors (SI, S2, and S3) consisted of two identical hybrids fixed back to back with 
perpendicular orientation of strips (to measure X and Y coordinates). Each silicon detector had an active volume of 
61 mm x 61 mm x 285 pm. The Si sensor was a single sided AC-coupled micro-strip detector. In the original design 
DU, strips had a pitch of 80 pm but since the required resolution for the PIENU experiment was of the order of 
300 pm, four strips were bound to one read-out line. For further reduction of the number of readout channels, the 
readout lines were interconnected with capacitors and only every fourth line was read out by an amplifier as shown in 
Figure[l0] Forty eight channels per silicon plane (288 channels in total) were read out. The capacitive network formed 
a charge division line where the reconstruction of the ionization amplitude and position was made by weighting of two 
or three channels that typically had a signal due to a charged particle traversing the sub-detector. The final position 
resolution (rms) provided by each silicon plane was 95 pm if at least two readout channels had a signal and 370 pm 
when the signal was recorded in only one readout channel. In order to reduce the data size, a hardware threshold was 
set to suppress channels with no hits. In SI and S2, the thresholds were optimized for the detection of pions while the 
S3 thresholds were set lower to ensure that the detection efficiency for decay positrons in at least one plane (X or Y) 
was higher than 99%. 
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Figure 10: Schematic drawing of the silicon strip detector readout scheme. 
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6.3. Performance 

This compact assembly of tracking devices was designed to minimize the major source of background for the 
measurement of 7T + —* e + v e decay due to decay of pions in flight. About 3.6% of the pions entering the PIENU 
detectoip~^|decay in flight upstream or inside the target. This type of decay, labeled pion-DIF for pion decay-in-flight, 
will have a low energy deposit in the target and can therefore mimic the energy deposition of a n + —> e + v e decay. 
About half of those pion-DIF events happen upstream of the target, between WC2 and SI, and can be identified by 
the tracking detectors as illustrated in Figure [TT] 



Figure 11: Schematic representation of pion-DIF and pion-decay at rest (PDAR) event. Events with decay in flight show an inconsistency between 
the directions of the track segments reconstructed in WC1 and WC2, and SI and S2, while track directions are consistent for decays at rest. 


The distributions of the angle (kink angle) between the track detected by WC1 and WC2 and the track detected 
by SI and S2 for different decay types are shown in Figure [12] for MC and on Figure [13] for data. Pion-DIF events 
between WC2 and SI can be distinguished due to their higher kink angle. 



Figure 12: Kink angle measured by WC1-WC2 and S1-S2 tracking 
devices from MC. About half of the pion-DIF events happen after SI 
and cannot be distinguished from pion-DAR (PDAR) or muon-DAR 
(MDAR) events. DAR stands for decay-at-rest. 



Figure 13: Kink angle for n + —» e + v e CIWa/( 7 Y) > 55 MeV) and n + —» 
H + —» e + (Tmj/( 7 Y) < 30 MeV) data events. The presence of pion-DIF 
events in the n + —» fi + —» e + sample can be clearly seen at angles > 
15°. 


14 This result is obtained from the Geant4 simulation. A pion is denoted as entering the detector if it enters the first detector element WC1 
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7. Trigger and data acquisition system 

7.1. n + —* e + v e run trigger 

Pions and a fraction of beam muons and positrons (used for the calibration of some detectors) were selected by the 
trigger using energy deposit information in B1. Definition of a beam particle was done by requiring the coincidence 
of the beam counters Bl, B2, and B3. A coincidence of T1 and T2 counters defined the decay-positron signal. A 
coincidence of beam pion and decay positron signals (PIE) within a time window of -300 ns to 500 ns was the basis 
of the main trigger logic. Since 7t + —> p 1 v fI decays happen much more often than n + — * e + v e decays, a Prescale trigger 
precisely selected only 1/16 of PIE events. Meanwhile, the n + —> e + v e events were enhanced by the Early and HE 
(High Energy) triggers. The Early trigger selected decays which happened between 2 ns and 40 ns (excluding prompt 
events) after the pion stop. Due to the 26 ns pion lifetime, more than 70% of the 7t + —> e + v e events happen within 
this time range. The HE trigger was based on a VME-module which summed the energy deposited in the calorimeters 
in real time. The trigger signal was produced by this module for events which had a high energy deposited in the 
Nal(Tf) and Csl spectrometers. The energy threshold was set 4 MeV below the highest energy of the positron from 
the n + —» p + — > e + decay chain. Almost all the n + —> e ' v,, events (with the exclusion of the tail events which extend 
below the e + energy spectrum from the n + —> /r 1 —> e + decay chain) were selected by this trigger. Those three triggers 
constituted the “physics triggers”. Additional triggers were used for data quality checks and calibration purposes. 
The Xe-lamp trigger provided flashes twice a second to all Csl crystals. Both these triggers ( Xe-lamp and Cosmic ) 
were intended for Csl calibration and monitoring. Finally, the Beam positron trigger triggered on one of every 32 
beam positrons to calibrate the Nal(Tf’) crystal. During a run, all 6 triggers were used and several of them could be 
satisfied at the same time. The rates of the various triggers are shown in Table [2] The trigger signal issued by any of 
the six triggers was then latched by the pion (V) and the positron (t e +) timings in order to reference the gate of the 
data acquisition modules to the incoming and decay particles, respectively. These latched signals triggered the data 
acquisition. T e + was used to trigger the data acquisition by the VME modules (VF48 and VT48, which are described 
below) while T-i triggered the COPPER data acquisition, described below. 


Table 2: Trigger rates 


Trigger 

Rate [Hz] 

Pion stop in target 

5xl0 4 

Physics triggers 

Early trigger 

160 

HE trigger 

170 

Prescale trigger 

240 

Other triggers 

Cosmics trigger 

15 

Beam Positron trigger 

5 

Xe lamp trigger 

2 

Total Triggers 

600 


7.2. Data acquisition system 

COPPER: This 500-MHz-sampling-frequency Flash ADC (FADC) system was based on the COPPER (COm- 
mon Pipelined Platform for Electronics Readout) platform. The detailed characteristics of the COPPER system have 
been described in Refs. lfl8l fl9l . The PIENU experiment was equipped with four COPPER board^Jto digitize the 
signals coming from the 23 PMTs reading out the plastic scintillators, analog sum of Nal(Tf) PMTs, and partial ana¬ 
log sums of Csl array PMTs. The time window of the signals recorded by COPPER covered 7.74 ps (1.35 ps after and 


15 An additional COPPER board was installed during the 2012 run to record partial analog sums of the Csl crystal PMTs. 
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6.39 //s before the pion stop) to be able to detect pre- and post-pile up particles. The waveforms recorded by COPPER 
were fitted, on an event-by-event basis, using templates obtained from a spline interpolation of the average PMT pulse 
shape. 
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Figure 14: Measured energy in Nal(TT) and Csl array for a portion of collected data. The vertical line denotes separation of collected events into 
low-energy (dominated by n + —> /i + v^ decays) and high-energy (dominated by n + —> e + v e decays and pile-up). 




e + time [ns] 


Figure 15: Time spectrum defined by the time difference between T1 
and B1 hits for low-energy events. Zero time corresponds to the time of 
the pion stop. The solid (red) line represents the fit function describing 
the 7r + —» fi + —» e + decay chain, the dashed (blue) line represents the 
sum of backgrounds, represented by two dotted lines which correspond 
to muon decays and pion decays in flight. 


Figure 16: Time spectrum defined by the time difference between T1 
and B1 hits for high-energy events. Zero time corresponds to the time 
of the pion stop. The solid (red) line represents the fit function describ¬ 
ing the n + —» e + v e time spectrum, the dashed (blue) line represents the 
n + —> n + —> e + decay chain events which due to pulse pile up ended up 
in the high energy spectrum, and the dashed line represents the sum of 
the remaining background mechanisms like charged pileup, radiative 
decays, etc. 
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VF48: The VF48 is a 48 channel 60 MHz FADC I l20l . It has a resolution of 10 bits and a dynamic range of 
+250 mV. All the Nal(Tf) and the Csl PMT signals as well as all the silicon strip channels were read out by VF48 
modules. This made a total of 404 channels (Nal(Tf): 19, Csl: 97, silicon strips: 288) read out by 10 VF48 modules). 
Typically, the VF48 read out 40 samples around the time defined by the decay positron. All VF48 modules were fed 
with the same 20-MHz clock provided by another VME module. This clock was multiplied internally in the VF48 to 
reach 60-MHz sampling. Data suppression was implemented in all channels except those reading the NalflY) crystal 
PMTs. Additionally, due to the slow NalfIT) response, its waveforms were digitized at a rate of 30 MHz instead of 
60 MHz. 

VT48: Discriminated signals from the three wire chambers (WC1, WC2, and WC3), all the logic signals issued by 
the PMT signals after discrimination, and most trigger logic signals were read out by multi-hit time-to-digital converter 
El. This device, based on the AMT-2 chips Il22ll . has a 25-MHz on-board clock which is multiplied to achieve 
0.625 ns resolution. All VT48s, however, were fed with an external 25-MHz clock to ensure the synchronization of 
modules. One board can read out 48 channels for up to 20 /-is. To optimize the dead-time, only two channels were 
read out for 20 fis to detect long-lifetime backgrounds originating from the beam while all other channels were read 
out 4.0 fjs before and after the time defined by the decay positron signal in Tl. In total, 11 VT48 modules were used 
in the experiment. 

The PIENU data acquisition system consisted of three VME crates (two VMEs hosted the VF48 and VT48 mod¬ 
ules while the third mostly ran slow control modules and modules used by the COPPER system) controlled by VME 
master modules, and four COPPER boards with a processor on each board. Each processor was running the asso¬ 
ciated front end programs to transfer the data to a host computer. Collection of the data was done by the MIDAS 
data acquisition system which incorporates an integrated slow control system with a fast online database and a history 
system ||23| . The MIDAS server computer can be controlled via a web interface. All the information and errors from 
the DAQ modules were issued on the web page and programs checking the quality of the data online were connected 
to MIDAS during data taking. 


8. Conclusion 


The PIENU detector was commissioned in 2009 and was operated during physics and supplemental measurements 
including NaKTf) response measurements until the end of 2012. The branching ratio R'^‘ will be determined from 
the observed numbers of n + —» e + v c and n + —> fi + —> e + events. The summed energy in the Nal(T) and Csl array 
for events in the time window -300 to +500 ns in relation to the pion stop time is shown in Fig. [Mjwhere the energy 
spectmm is divided into two regions at 52 MeV: low-energy events dominated by n + —> /C —> e + decays, and high- 
energy events dominated by n + —> e + v,, decays. The time spectra for these two regions shown in Figs. 15 and 16 
are used to determine the raw branching ratio taking into account various background components such aspile-up hits 
from previous muon decays, and pion decays-in-flight. The R elfJ is then determined after various corrections are made 
(See Section 1). The excellent performance of the crystal array and tracking array, paired with good time resolution 
achieved in plastic scintillators (Bl, Tl) were important aspects of the PIENU apparatus ensuring the possibility of a 
successful measurement of R e ^‘ and improved searches for heavy neutrinos 0. 
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